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Abstract 

The low-frequency optical excitations of AA-stacked bilayer graphene are investi- 
gated by the tight-binding model. Two groups of asymmetric LLs lead to two kinds 
of absorption peaks resulting from only intragroup excitations. Each absorption peak 
obeys a single selection rule similar to that of monolayer graphene. The excitation 
channel of each peak is changed as the field strength approaches a critical strength. 
This alteration of the excitation channel is strongly related to the setting of the Fermi 
level. The peculiar optical properties can be attributed to the characteristics of the 
LL wave functions of the two LL groups. A detailed comparison of optical properties 
between AA-stacked and AB-stacked bilayer graphenes is also offered. The compared 
results demonstrate that the optical properties are strongly dominated by the stacking 
symmetry. These optical properties could be verified by optical measurements. 
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1 Introduction 

Few-layer graphenes (FLGs) are very exotic nanomaterials owing to their nanoscale inter- 
layer distance, hexagonal symmetry, and stacking configurations. FLGs have attracted nu- 
merous investigations on band structures, optical spectra, ^'^"^^ electronic excitations, ^'^^'^ 
and transport properties. The presented properties seem to make FLG-based materials 
excellent candidates for application in electronic and photonic devices. Since the phys- 
ical properties of FLGs are strongly affected by the stacking configurations, FLGs with 
different stacking configurations have attracted considerable experimental and theoretical 
research. Monolayer (MG), AA-stacked bilayer (AABG) and AB-stacked bilayer (BEG; bi- 
layer Bernal) graphenes are three prototypical FLGs. For the low-lying energy dispersions, 
MG exhibits isotropic hnear bands near the Fermi level {Ep — 0); these bands become 
gradually anisotropic parabohc bands as the energy exceeds the region of ±0.5 eV.^ In 
AABG, the linear bands in MG change into two pairs of linear subbands with slightly 
different slopes.^ These two pairs cross at Ep — and are asymmetric about Ep — 0. For 
BBG, two pairs of parabolic subbands exist that are asymmetric about Ep — 0.^'^^ The 
conduction and valence bands of the first pair slightly overlap about Ep — 0. 

In the presence of a uniform perpendicular magnetic field B = Bqz, the zero-field energy 
bands of FLGs become the dispersionless Landau levels (LLs).^"'''^"'^ In MG, the low- lying 
LLs are characterized by the special relation E^'" oc y/n^''"Bo, where n'^ {n") is the quantum 
number of the conduction (valence) LLs. This special relation is broken as the energy 
exceeds the region of ±0.5 eV since the hnear dispersions become gradually parabolic 
as the energy exceeds this region.^ For AABG and BBG, two groups of asymmetric LLs 
exist. ^'^"^ These LLs are distributed away from a certain energy value in each group. The 
LL energies of both bilayer graphenes do not exhibit a simple relation similar to that of 
MG. These main features of FLGs would be refiected in the magneto-optical absorption 



3 



spectra. MG and BBG display different optical properties, e.g., different field-dependent 
absorption frequencies and distinct optical selection rules. ^'^^ 

In this work, the magneto-optical absorption spectra of AABG are calculated by gra- 
dient approximation^'-'^^'-'^^'^^ within the tight-binding model (TB). Two groups of LLs are 
divided based on the characteristics of the LL wave functions. The LL wave functions are 
clearly depicted and utilized to explain the main features of the optical absorption spectra. 
The results show that two kinds of absorption peaks exist in the absorption spectra. Each 
peak obeys a single selection rule. The excitation channel associated with each peak varies 
with a changing field strength. The optical properties can be reasonably comprehended by 
the LL spectra and the characteristics of the LL wave functions. A detailed comparison 
between AABG and BBG reveals that they possess different magneto-optical absorption 
spectra, reflecting the influences of stacking conflgurations on the electronic properties. In 
other words, this method offers another way to distinguish AABG from BBG in addition 
to the STM images. Moreover, it may also be used to discriminate AABG from MG, which 
can be hardly done by STM.^^ 



2 Landau level spectrum 

The geometric configuration of AABG is shown in Fig. 1(a). The primitive unit cell 
consists of four sublattices, Ai, Bi, A2, and B2. The subscripts 1 and 2 are, respectively, 
the indices of the first and second layer. Three atomic hopping integrals, ckq (=2.569 eV), 
Q!i (=0.361 eV), and (=-0.032 eV), are taken into account in this work. The first integral 
is the nearest-neighbor hopping integral on the same layer and the second and third are the 
interlayer interactions, as indicated in Fig. 1(a). B induces a periodic Peierls phase related 
to the vector potential A(r) = {0,Bqx,0). The Hamiltonian under a magnetic field is 
Hb — (P — eA(r)/c)^/2m-|-y(r), where m, P (= ^k), and V{r) are the electron mass, the 
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crystal momentum, and the lattice potential, respectively. Under the periodic condition, 
the primitive unit cell is enlarged^'^'^ and composed of four effective sublattices (denoted 
Ai, Bi, A2, and B2 for convenience, similar to the symbols of the zero- field wave functions) 
including 2Rb Ai, 2Rb Bi, 2Rb A^, and 2Rb B^ atoms, respectively. Rb is defined by 
Rb = $o/(^-Bofc^) = ^^^^p, where $0 is the flux quantum and b is the C-C bond length. 
Rb is inversely proportional to Bq and related to the dimension of H^;^'^'^ for example, Rb 
is 1975 for Bq — 40 T. That is to say, each LL wave function is the hnear combination of the 
four magnetic TB functions associated with the four effective sublattices. The Hamiltonian 
matrix elements in the presence of a magnetic field are 



where Ri,M is the position vector of A^^m or B^^m- 7s's(Ri,M, R-i'j') indicate the atomic 
interactions between the atoms at Ri,M and Ri',M', i-e., they are ao, ai, and in this 
work. The representation Ai^M (Bi^m) indicates the Mth (M = 1, 2...2Rb) A (B) atom on 
the ith (i = 1, 2) layer. 

The first Brillouin zone of AABG is shown in Fig. 1(b). For the low-energy electronic 
structure of AABG, the linear subbands of MG change into two pairs of linear subbands 
owing to the AA-stacking configuration and interlayer interactions, as shown in Fig. 1(c). 
The analytical solution of the energy dispersions of the first (second) pair can be described 



as E"''' = -ai ± I ( ao - ^3) bk + ^ (E^'" = ai ± | (ao + ^a) bk + ^),^ where E" and 



E'" are respectively associated with the conduction and valence bands and are crossing at 
-ai + ™ ~ -0.366 eV (cvi + ^ ~ 0.357 eV), as indicated by the black (red) lines, k is 
the wave vector measured from the K point. The conduction and valence bands of each pair 
are symmetric about the crossing energy. Furthermore, the two pairs of subbands intersect 
at Ep = and their slopes are slightly different. This causes the occupied and unoccupied 




R^,M)]dX], (1) 
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states to be asymmetric about = 0. For the wave functions of zero-field subbands (not 
shown), the first and second pairs show the special relations, A1/A2 — B1/B2 — —1 and 
A1/A2 — B1/B2 — 1, respectively. These relationships mainly result from the nearest- 
neighbor interlayer interaction ai, which leads to the separation of the first pair from the 
second pair. These main features of the zero-field subbands would be reflected in a LL 
spectrum. 

The magnetic field quantizes the two pairs of zero-field linear bands into fourfold de- 
generate LLs, as shown in Fig. 2(a) for 40 T. Based on the characteristics of the wave 
functions (Fig. 2(b)) discussed below, these LLs can be further divided into two groups of 
LLs, ^**LLs and ^"'^LLs. The first (second) group corresponds to the first (second) pair of 
zero-field subbands and is distributed away from the onset energy -0.366 eV (0.357 eV), as 
shown in Fig. 2(a) by the black (red) lines. The onset energy of the first (second) group 
is located at the crossing energy of the first (second) pair of zero-field subbands. The LLs 
of each group are symmetric about the onset energy of this group. However, the occupied 
and unoccupied LLs are asymmetric about Ep — 0. The LL energies of the two groups can 
be written as 

^ -ai + ^ ± ^ ( ao - 0^3) hJ 2eB,nr''" Iti (2a) 
£; ^ ai + ^ ± ^ ( «o + «3) hJ2eBonr^'f^/h, (2b) 

where rf{"'^^^ and 712"'^^^ are effective quantum numbers of the LLs in the first and second 
groups and defined below. The low-lying LL energies in each group linearly depend on \fB 
similar to the relationship of MG.^'^ 

The four degenerate wave functions of a LL are similar and thus only one of them is 
discussed here. Each LL wave function of AABG is linearly combined by four TB functions 
associated with the four effective sublattices. These TB functions display oscillatory modes 
and localizaed features. The four TB functions localized at a certain position are chosen 
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for discussions (Fig. 2(b)). ^'^"^ Through appropriate fitting, the wave functions of the nth 
^^%L (^""^LL) can be expressed as Ai — —A2 oc (pn-i (x) and Bi — —B2 oc ipn-2 (x) 
{Ai — A2 (X (fin-i {x) and Bi — B2 oc (pn-2 {x))- <Pn {x) is the product of the nth-order 
Hermite polynomial and Gaussian function, ^'^'^^'^"^ where n is the number of zero points of 
(fin (x) and chosen to define the quantum number of a LL.^^-^ For convenience, the zero-point 
numbers of the A atoms are chosen as effective quantum numbers of the LLs. n\^^^s and 
rfj^''s are, respectively, the effective quantum numbers of the LLs in the first and second 
groups, as shown in Fig. 2(a). Thus, the effective quantum numbers of the nth conduction 
(valence) LLs in both groups are n — 1, which are similar to those of the LL wave functions 
of MG.^ Furthermore, the LL wave functions in Fig. 2(b) show a special relationship 
similar to those of the zero-field wave functions, i.e., A1/A2 = B1/B2 = —1 and A1/A2 — 
B1/B2 — 1 for the first and second groups, respectively. These relations mainly originate 
from the specific stacking configuration and should strongly affect the optical properties, 
e.g., the optical- absorption peak structure and the optical selection rules. The LL wave 
functions in both groups exhibit features similar to those of MG, i.e., they display similar 
oscillation modes, localization features, and a combination of the A and B atoms in each 
layer. Obviously, the characteristics of the LLs reflect those of the zero-field subbands, 
i.e., the existence of the two LL groups, the specific onset energies of the two groups, the 
symmetric structure about the onset energy of each group, the asymmetry of the occupied 
and unoccupied LLs about Ep — 0, and the specific relations of the LL wave functions of 
the two groups. 



3 Magneto-optical properties 

The main features of the LL spectra would be reflected in optical excitations. At zero 
temperature, there exist only excitations from the occupied to the unoccupied states. Based 
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on Fermi's golden rule, the optical absorption function is given by 
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xlm{ \ , \\ (3) 



Eh' (k, nl'f"^ - Eh (k, n1$"^ - u - iT 
E is the unit vector of an electric polarization and E || x is taken into account in this 
work, r (=1 meV) is a broadening parameter and often affected by temperature and defect 
effects, h and h' represent the occupied and unoccupied states, respectively. {'^^' {n^^2^ 1^)1^- 
P/me|'?/'^(nJ^2; k)) is the velocity matrix element (denoted M^'^'j derived from the dipole 
transition and calculated by gradient approximation. ^^'^^ Through detailed calculations, 
Mh'h is expressed as 

J2 E [Al'* xB^]V^{A,My.\HB\Bj,M'y.)+h.c. (4) 

i,j=l,2 M,M'=l 

Eq. (4) corresponds to the two hopping integrals ccq (the terms for i — j) and 0:3 (the terms 
for i 7^ j), and the terms associated with aQ dominate the value of Eq. (4). For the sake 
of convenience, the dominant term Af'* x is represented by Mj^''^(q;o) in the following 
discussions. 

The low-frequency optical absorption spectra for Bq — 0, 20 T, and 40 T, are shown in 
Fig. 3(a). The zero- field spectrum does not show any peak structure and the absorption 
rate is zero as the energy is below 0.71 eV. The vanishing absorption rate corresponding 
to the intersubband excitations is owing to the special relations of the wave functions in 
the two subbands. The emergence frequencies of the intrasubband excitations associated 
with the first and second pairs are nearly 0.73 eV (~ 2ai — 2^^; indicated by the black 
dot) and 0.71 eV (~ 2ai + 2^; indicated by the red dot), respectively. For = 20 T 
and 40 T, absorption peaks basically result from the excitations between two LLs in the 
same group, i.e., only intragroup excitations exist in the absorption spectrum. Each peak 
can clearly be identified. The excitation channels of the first (second) kind of absorption 
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peaks, nl'"'^^^ n^'"'^-''^ {nl'"'^^^ ria''''''^-'^), are indicated by black (red) dots in Fig. 3(a). 
^c,v,eff _^ ^/c.'y.e// ^^c,v,eff _^ ^ic,v,eff^^ represents the excitation channel from the occupied 

i^*LL with nl'"'^^^ (^"'^LL with n/'^") to the unoccupied ^^*LL with n'^'"'^^^ (2"<^LL with 
fi'^'^'^ffy The mth absorption peak frequency of the first (second) kind is denoted u}"{ 
(u;^). For = 20 T (40 T), ul^ and 00^2 originate from 7" 8^ and 8^ T (3^ A" and 
4" 3*^) respectively, i.e., they originate in the excitations from the occupied conduction to 
the unoccupied conduction LLs and from the occupied valence to the unoccupied conduction 
LLs, respectively. uj\'y and ijj\2 are merged owing to their almost identical frequencies and 
strongly affected by the setting of Fermi energy (the Fermi energy is set to Ep = 0). Except 
for cull and 00^2, the other peaks, ou^ and ou^ for m > 2, come from the excitations between 
the occupied valence and unoccupied conduction LLs. For Bq = 20 (40) T, the excitations 
7" 8" and 8^ 7" (3^ A" and 4^ 3^) result in uj^ and respectively. For 20 T, the 
other channels, (m + 5)^ — >■ (m + 6)*^ and (m + 6)^ — >■ (m + 5)*^ possess the same frequency 
in the first (second) group which result in cu^ ((^S)- However, cu'^ (f^S) for 40 T originates 
from (m + 1)^ — >■ (m + 2y and (m + 2)^ — >■ (m + 1)'^ possessing the same frequency in the 
first (second) group. Obviously, the excitation channels of the mth peaks related to the 
different field strengths may be different, cj^ and o;^ are also slightly different. The former 
is higher than the latter, and they form a pair-like structure. This reflects the asymmetric 
structure of the LLs. The vanishing peak structure in the low-frequency region below a;|2, 
except for ouli and UJ22, can be ascribed to the characteristics of the LL wave functions 
and is discussed below. Simply said, the magneto-optical absorption spectra reflect the 
main features of the zero-fleld spectrum, i.e., two kinds of absorption peaks exist and only 
intragroup excitations are allowed. 

The optical selection rules for cun and UJ22 can be represented by | Ann | ( 
I An22| (= 712'"^^ - nY^^ )= 1- These rules are similar to those of the LLs in MG. The 



h',ejf h,eff 

n-i — n-i 
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selection rules can be comprehended by the characteristics of the LL wave functions. In Eq. 
(3), Mj^''*(q;o) dominates the excitations of the absorption peaks. M!^''^{ao) (= Af'* x B^) 
has non-zero values only when Ai and Bi {A2 and B2) own the same (pn{x) because of 
the orthogonality of (pn{x)- Since Ai {A2) in the nth LL and Bi {B2) in the n + 1th LL 
for both groups own a same (fin{x), the selection rules |Anii| = |An22| = 1 can be easily 
obtained. Furthermore, the disappearance of the intergroup excitations can be ascribed to 
the fact that the two products, Ai x Bi and A2 x B2, in the intergroup excitations cancel 
each other out due to the special relationship of the wave functions. 

The field-dependent absorption frequencies associated with the first and other mth 
peaks are shown in Figs. 3(b) and 3(c), respectively. In Fig. 3(b), at Bq — 60 T, the 
excitation channel of (xil^ (1^22) comes from 2*^ 3^^ (3" 2"). With decreasing field 
strength, the frequencies of ujI^ and UJI2 approach each other and merge at a sufficiently 
small field strength. Furthermore, the excitation channels of the first peak are altered as 
the field strength decreases to a critical field strength. For example, the channel 2'^ 3'^ 
(3" 2") changes into 3" ^ 4^ (4^ ^ 3^^) and 4^ ^ 5" (5^ ^ 4^) at = 47 T and 35 
T (indicated by two yellow lines), respectively. As Bq is reduced, the excitation channels 
become those associated with the LLs possessing larger effective quantum numbers. The 
main reason for this is that the setting of the Fermi level strongly affects which LL is con- 
sidered the highest occupied (lowest unoccupied) one. The discontinuity of field-dependent 
absorption frequencies can be used in optical experiments to determine the excitation chan- 
nels. Moreover, the excitation frequency of each excitation channel in its existent region 
is proportional to i?o, a behavior similar to the excitation frequencies of MG. The other 
excitation channels shown in Fig. 3(c) show features similar to those in Fig. 3(b). The 
convergent frequencies of cun and a;22 at the weak field strength are approximately 0.73 eV 
and 0.71 eV and correspond to the two emergence frequencies in the zero-field absorption 
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spectrum, respectively. Since the absorption peak intensities related to the LLs are strong, 
optical measurements can reasonably determine the values of ai and through observing 
the convergent frequencies of absorption peaks. 

In addition to AABG, the magneto-optical properties of BBG are also discussed in a 
previously pubhshed work.^^ AABG and BBG show similar LL spectra, e.g., two groups 
of LLs and an asymmetric structure. However, these two prototypical bilayer graphenes 
display totally different optical properties. AABG exhibits two kinds of absorption peaks, 
and only intragroup excitations that follow a single optical selection rule take place. BBG 
produces four kinds of absorption peaks, and both intra- and inter-group excitations that 
follow complex optical selection rules take place. These optical properties can be compre- 
hended by obtaining the characteristics of the LL wave functions. The differences between 
the optical absorption spectra of AABG and BBG imply that the optical properties can 
reflect the influences of different stacking configurations. The dissimilarities between these 
two graphenes are helpful for distinguishing AABG from BBG via optical measurements. 
Furthermore, both the magneto-optical properties of AABG and BBG reflect the main 
features of zero-field optical properties. 

4 Conclusions 

In summary, the TB calculations show that the LLs are asymmetric about Ep — and can 
be divided into two groups based on the characteristics of the wave functions. These two 
groups lead to two kinds of optical-absorption peaks Wn's and a;22's associated with only 
intragroup excitations of the first and second groups, respectively. The absorption frequen- 
cies of cu^^s and a;^'s are slightly different and form pair- like structures, which originate 
from the asymmetry of LLs. The optical selection rules can be reasonably explained by the 
characteristics of the LL wave functions. The selection rules of the two kinds of peaks are 
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Anil = An22 = ±1 similar to that of MG. The similar selection rules of AABG and MG 
mainly originate from the resembling characteristics of their LL wave functions. Further- 
more, each a;^'s (o'^'s) corresponds to different excitation channels within the different 
field strength region. The field-dependent absorption frequencies for each excitation chan- 
nel are hnearly dependent on \/B^, and resemble the absorption frequencies of MG. The 
convergent absorption frequencies at the weak field strength might be helpful and reliable 
in determining the interlayer atomic interactions cti and a^. The different magneto-optical 
properties of AABG and BEG reflect the influences of different stacking configurations. 
The differences between the two prototypical bilayer graphenes can help experimental re- 
searchers discriminate AABG from BBG. The above-mentioned magneto-optical properties 
could be confirmed by magneto- absorption spectroscopy measurements. 
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5 Appendix 



In the absence of external fields, the Hamiltonian is Hq — P^/2m + V{r), where m, P 
(= hk), and V{r) are the electron mass, the crystal momentum, and the lattice potential, 
respectively. In the presence of a uniform magnetic field B = Bqz, the Hamiltonian is 
i?B — (P — eA(r)/c)^/2m + V{r), where A(r) = (0, SqX, 0) is the vector potential. A 
Peierls phase induced by the magnetic field leads to a periodic condition and thus cause 
the primitive uint cell to be enlarged, i.e., the primitive unit cell in the absence of external 
fields is composed of four atoms, Ai, Bi, A2 and B2, while the enlarged primitive unit 
cell in the presence of B comprises 2Rb Ai, 2Rb -Bi, 2Rb A2, and 2Rb B2 atoms. Rb = 
^o/i^^Bob'^) = ^ is associated with the dimention of the Hamiltonian matrix. The 
Peierls is defined by 

AG^ [ {Ri^j, -Rij) ■ A[Ri,j, + X{Ri,,f -Rij)]dX, (A.l) 
Jo 

where Rjj is the position vector of Aij or Bij for i = 1,2 and j — 1,2, ...2Rb. In the 

sequence of the bases: \A\), \Al), \Bl), \Bl); \Al), \Al), \Bl), \A\j^^), \Al^^), 

\B2Rb)i \^2Rb)i Hamiltonian matrix related to the enlarged unit cell is expressed by 



H 



Hs,j=i 


Ha,j=2 
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Ha,j=3 














Hl,j=3 


Hs,j=3 


Ha,j=4 
















Hs,j=A 




Ha,j=2RB-l 















Hij=2RB-l 


Hs,j=2RB-l 


Ha,j=2RB 
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^^a,j=2RB 


Hs,j=2RB 



\ 



(A.2) 
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where 



Ha^a^ Ha^.a^. HaIb^ HaIb^ 

3 3 3 3 3 3 3 3 



\ 
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J "J 



if^2^i Ha^a^ Ha^b^ Ha^b- 

3 3 3 3 3 3 ' ■ 
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3 3 3 3 3 3 
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3 3 3 3 
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3 3 
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/3i 




, (A.3) 



and 
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y y 



iij and q are expressed as 



iij = exp{^[-/c^- - + TT— (j - 1 + -)]} 

+ exp{^[-A;^- + /c^^ - tt— (j - 1 + -)]}, 
g = exp[i/ca;6]. 



(A.4) 



(A.5) 



(A. 6(a)) 
(A. 6(b)) 



The optical absorption rate is dominated by the velocity matrix M^'^. Through the 
gradient approximation, the velocity matrix is simplified as the product of three matrices, 
the initial state (occupied state; V^'*), final state (unoccupied state; V''*'), and the first-order 
differential of the Hamiltonian matrix element versus the wave vector k (Vk-f^ij)- The last 
term corresponds to the direction of electric polarization, i.e., it is dHij/dk^ {dHij / dky) 
for polarization along the x-direction (y-direction) . The elements in the third matrix are 
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non-zero only for H^j associated with the hopping integrals, as are the velocity matrix 
elements. In other words, when the velocity matrix element does not vanish, the initial 
and final states in the product should be the two states in a hopping process corresponding 
to non-vanishing hopping integrals. In this work, the polarization is along the x-direction 
and M^'^ is expressed as 



2Rb 



d/dK [i^i'Mk I Hb I ^'Im'i 

M,M'=1 



^P'l)+h.c., (A. 7) 



where ip is the wave functions of the A or B atom. Eq. (A. 7) is associated with the three 
hopping integrals ao, cti and 0:3. The term related to ai is vanishing since the relative 
differential value is zero. M'^''^ can be simplified as 



1=1,2 



2Rb 



\Hb\B 



M=l 



) + h.c. 



i=l,2; i^j 



2Rb 



J2 dH/dK{A i,Mk \Hb \ Bj^M'k) 
M,M'=l 



B^)6M,M'±i + h.c. (A. 8(a)) 



(A. 8(b)) 



The first term, M^^^{aQ), is associated with ckq and the second term, Ml^'^{as), corresponds 
Qis; Mj^''*(q;o) dominates the value of the velocity matrix. 
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Figure Captions 



FIG. 1. (a) The geometric structure, (b) the first Brillouin zone, and (c) the low-lying sub- 
bands of AA-stacked graphene. ckq (2.598 eV) is the nearest-neighbor hopping integral 
and the two important interlayer interactions are cti (0.361 eV) and (—0.032 eV). 

FIG. 2. (a) Landau levels and (b) Landau level wave functions of AA-stacked graphene at 
Bq = 40 T. n^'^^'s and n2^^s are the effective quantum numbers of the first and second 
group of Landau levels, respectively. 

FIG. 3. (a) The optical absorption spectra of AA-stacked graphene at Bq = 0, 20, and 40 T. 
The field-dependent absorption frequencies associated with (a) the first and (b) other 
mth absorption peaks. The symbols n ^ n' describe the excitation channels of the 
absorption peaks. 
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